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Reduction operators (called also nonclassical or Q-conditional symmetries) 
of variable coefficient semilinear reaction-diffusion equations with exponential 
source f(x)ut — (g(x)u x ) x + h(x)e mu are investigated using the algorithm in- 
volving a mapping between classes of differential equations, which is generated 
by a family of point transformations. A special attention is paid for checking 
whether reduction operators are inequivalent to Lie symmetry operators. The 
, derived reduction operators are applied to construction of exact solutions. 
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1 Introduction 



Various processes in nature are successfully modeled by nonlinear systems of par- 
tial differential equations (PDEs). In order to study behavior of these processes, 
it is important to know exact solutions of corresponding model equations. Lie 
symmetries and the classical reduction method present a powerful and algorith- 
mic technique for the construction of exact solutions (of systems) of PDEs |14pi6j . 
In [2] Bluman and Cole introduced a new method for finding exact solutions of 
PDEs. It was called "non-classical" to emphasize its difference from the classical 
Lie reduction method. A precise and rigorous definition of nonclassical invari- 
ance was firstly formulated in [6] as "a generalization of the Lie definition of 
invariance" (see also |26j). Later, operators satisfying the nonclassical invariance 
criterion were called, by different authors, nonclassical symmetries, conditional 
symmetries and Q-conditional symmetries [2J|31|7|[T2]. Until now all names are 
in use. See, e.g., 1 1 1|. fT5 . 20 for comprehensive reviews of the subject. Follow- 
ing Ref. [19] we call nonclassical symmetries reduction operators. The necessary 
definitions, including ones of equivalence of reduction operators, and statements 
relevant for this paper are collected in 
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The problem of finding reduction operators for a PDE reduces to integration of 
an overdetermined system of nonlinear PDEs. The complexity increases in times 
in the case of classification problem of reduction operators for a class of PDEs 
having nonconstant arbitrary elements. 

The experience of classification of Lie symmetries for classes of variable coeffi- 
cient PDEs shows that the usage of equivalence and gauging transformations can 
essentially simplify the group classification problem and even be a crucial point 
in solving the problem [8, , 221124] . The above transformations are of major impor- 
tance for studying reduction operators since under their classification one needs to 
overcome much more essential obstacles then those arising under the classification 
of Lie symmetries. 

In [23] we propose an algorithm involving mapping between classes for finding 
reduction operators of the variable coefficient reaction-diffusion equations with 
power nonlinear ity 



where /, g and h are arbitrary smooth functions of the variable x and m is an 
arbitrary constant such that fgh ^ and m/0,1. In [23] reduction operators of 
the equations from class (pQ) with m ^ 2 were investigated using this algorithm. 
The case m = 2 was not systematically considered since it is singular from the Lie 
symmetry point of view and needs an additional mapping between classes (see [21] 
for more details). Nevertheless, all the reduction operators constructed in [23] for 
the general case m/0,1,2 are also fit for the values m = 0, 1, 2. 

In this paper we implement the same technique to find reduction operators of 
the variable coefficient reaction-diffusion equations with exponential nonlinearity 



Here /, g and h are arbitrary smooth functions of the variable x, fgh ^ and m 
is an arbitrary nonvanishing constant. 

The structure of this paper is as follows. For convenience of readers section [2] 
contains a short review of results obtained in [21J and used here. Namely, in 
this section the necessary information concerning equivalence transformations, the 
mapping of class @ to the so-called "imaged" class and the group classification of 
equations from the imaged class is collected. Moreover, all additional equivalence 
transformations connecting the cases of Lie symmetry extensions (cf. table 1) are 
first found and presented therein. As a result, the classifications of Lie symmetry 
extensions up to all admissible point transformations in the imaged and, therefore, 
initial classes are also obtained. Section [3] is devoted to the description of the 
original algorithm for finding reduction operators of equations from class (j2]) using 
a mapping between classes generated by a family of point transformations. The 
results of section [4] are completely new and concern the investigation of reduction 
operators for equations from the imaged class. Reduction operators obtained in 
an explicit form is used for the construction of exact solutions of equations from 
both the imaged and initial classes. 




(1) 




(2) 
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2 Lie symmetries and equivalence transformations 

Class (|2]) has complicated transformational properties. An indicator of this is 
that it possesses the nontrivial generalized extended equivalence group, which 
does not coincide with its usual equivalence group, cf. Theorem Q] below. To 
produce group classification of class ([2]), it is necessary to gauge arbitrary elements 
of this class with equivalence transformations and a subsequent mapping of it 
onto a simpler class j2TJ[2l]- It appears that the preimage set of each equation 
from the imaged class is a two-parametric family of equations from the initial 
class ([2]). Moreover, preimages of the same equation belong to the same orbit of 
the equivalence group of the initial class. It allows one to look only for simplest 
representative of the preimage to obtain its symmetries, exact solution, etc., and 
then reproduce these results for two-parametric family of equations from the initial 
class using equivalence transformations. 

Theorem 1. The generalized extended equivalence group G~ xp of class $2$) consists 
of the transformations 

i=5it + 52, x = (p(x), u = 63U + ifrix), 

f= J, 9 = o (p x g, h = e 5 3 h, m = —, 

where <p is an arbitrary non- constant smooth function of x, if) = 64 J + £5, 
and Sj, j = 0, 1, . . . , 5, are arbitrary constants such that (5o<5i^3 7^ 0. 

Corollary 1. The usual equivalence group of class ([2]) is the subgroup of G*^ p 
singled out by the condition 64 = 0. 

The transformations from associated with varying the parameter 5q in fact 
do not change equations from class ([2]) and hence form the gauge equivalence group 
of this class. The values of arbitrary elements connected by a such transformation 
correspond to different representations of the same equation. 

Analogously to the power case we at first map class ([2]) onto its subclass 

f(x)u t = (f(x)u x ) x + h(x)e u (3) 

(we omit tildes over the variables) , using the family of equivalence transformations 
parameterized by the arbitrary elements /, g and m, 



t = sign(/(x)g((x))t, x 



dx, u = mu. (4) 



The new arbitrary elements are expressed via the old ones in the following way: 

i 

~ 1 ~ g(yX^j 

fix) =g(x) = sign(g(x))\f(x)g(x)\i, h{x) = m — -r 

The next step is to change the dependent variable in class ([3]): 



hix), m = 1. 



v(t,x) = u(t,x) + G(x), where G(x) = In \f(x) l h{x)\. (5) 
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Finally we obtain the class 

v t = v xx + F(x)v x +se v + H(x), (6) 

where e = sign(/(x)/i(x)) and the new arbitrary elements F and H are expressed 
via the arbitrary elements of class ([6]) according to the formulas 

F = fxf~ 1 , H = —G xx — G X F. (7) 

All results on Lie symmetries and exact solutions of class © can be extended 
to class ([3]) by the inversion of transformation ([5]). 

The arbitrary elements / and h of class ([6]) are expressed via the functions F 
and H in the following way: 

/ = Co exp (J" Fdx) , h = eco exp (J" Fdx + G) , 
where G = j e ~f Fdx (c\ — j He$ Fdx dx S j dx + C2- 

Here Co, c\ and C2 are arbitrary constants, cq ^ 0. The constant Co is inessen- 
tial and can be set to the unity by an obvious gauge equivalence transformation. 
The equations from class ([3]), which have the same image in class ([6]) with re- 
spect to transformation ©, i.e., whose arbitrary elements are given by © and 
differ only by values of constants c% and C2, are G^p-equivalent. The equivalence 
transformation 

t = t, x = x, u = u + c\ j e~ f Fdx dx + C2 (9) 

maps an equation ([6]) having / and h of the form ([8]) with cf + 7^ to the 
one with c\ = C2 = 0. Hence up to G^ p -equivalence we can consider, without 
loss of generality, only equations from class ([3]) that have the arbitrary elements 
determined by (JSJ) with c\ = C2 = 0. 

Theorem 2. The generalized extended equivalence group G~ of class ([6]) coicides 
with its usual equivalence group and is formed by the transformations 

i = S^t A- 82, x = Six + 63, v = v — \n.5l, 
F = 5{ l F, H = &{ 2 H, 

where 6j, j = 1,2,3, are arbitrary constants, 5\ 7^ 0. 

The kernel of the maximal Lie invariance algebras of equations from class ([6]) 
is the one-dimensional algebra (c\). It means that any equation from class ([6]) is 
invariant with respect to translations by t, and there are no more common Lie 
symmetries. 

Theorem 3. G~ xp -inequivalent cases of extension of the maximal Lie invariance 
algebras in class ([6]) are exhausted by ones presented in table d 
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Table 1. The group classification of the class vt = v xx + F(x)v x + ee v + H(x). 



N 


F{x) 


H(x) 


Basis of A max 





V 


V 


St 


1 


ax~ 1 + fix 


j3x~ 2 + 2/i 


<9 t , e- 2 ^(d t - fj,xd x + 2fid v ) 


2 


ax^ 1 




d t , 2td t + xd x - 2d v 


3 


fix 


7 


d u e~^d x 


4 


X 


7 


dt, d x 


5 


fix 


2/i 


d t , e~^d x , e- 2 ^(d t - fixd x + 2fxd v ) 


6 


X 





d t , d x , 2td t + (x - Xt)d x - 2d v 



Here A € {0, 1} mod G~ p , \i = il mod G~ p ; a,/3,7 are arbitrary constants, 
a 2 + (3 2 / 0. We also have 7 7^ 2\i and 7 7^ in cases 3 and 4, respectively. 

The corresponding results on group classification of class ([2]) up to G^-p-equiv- 
alence is given in table 3 of |21j . 

Additional equivalence transformations between G~ -inequivalent cases of Lie 
symmetry extension are also constructed. The pairs of point-equivalent cases from 
table Q] and the corresponding transformations are exhausted by the following: 

Ih+2, 5^6k_ n : t = — e 2 ^, x = e> lt x, v = v-2iit, 

( 10 ) 

4 1— >• 4|^_ , 61— >-6|^ =Q : t = t, x = x + Xt, v = v. 

The inequivalence of other different cases of table 1 can be proved using differences 
in properties of the corresponding maximal Lie invariance algebras, which should 
coincide for similar equations. Thus, the dimensions of the maximal Lie invariance 
algebras are one, three and two in the general case (case 0), cases 5 and 6 and 
the other cases, respectively. In contrast to cases 1-3, the algebra of case 4 is 
commutative. The derivative of the algebra of case 3 has the zero projection on 
the space of t, and this is not the case for cases 1 and 2. Possessing the zero 
(resp. non-zero) projection on the space of t is an invariant characteristic of Lie 
algebras of vector fields in the space of the variables t, x and v with respect to 
point transformations connecting a pair of evolution equations since for any such 
transformation the expression of the transformed t is well known to depend only 
on t PUS]. 

More difficult problem is to prove that there are no more additional equiv- 
alences within a parameterized case of table 1. (In fact, all the cases are pa- 
rameterized.) This needs at least a preliminary study of form-preserving |9j (or 
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admissible |17yi8j) transformations. In contrast to transformations from the cor- 
responding equivalence group, which transform each equation from the class C of 
differential equations under consideration to an equation from the same class, a 
form-preserving transformation should transform at least a single equation from C 
to an equation from the same class. The notion of admissible transformations is a 
formalization of the notion of form-preserving transformations. The set of admis- 
sible transformations of the class C is formed by the triples each of which consists 
of the tuples of arbitrary elements corresponding to the initial and target equa- 
tions and a point transformation connecting these equations. It is obvious that 
each transformation from the equivalence group generates a family of admissible 
transformations parameterized by arbitrary elements of the class C 

A preliminary description of the set of admissible transformations of the class © 
which is sufficient for our purpose, is given by the following statements. 

Proposition 1. Any admissible point transformation in the class ([6]) has the form 

t = T(t), x = 5y / T t x + X(t), v = v-lnT t , 

where 5 = ±1 and T and X are arbitrary smooth functions oft such that T > 0. 
The corresponding values of the arbitrary elements are related via the formulas 

s/T t Tt 1 T t T t 2 

Corollary 2. Only equations from the class © whose arbitrary elements have 
the form 

F = /UX + A +^-, = 7+ 7 f , 2 , (11) 

X + K \X + KY 

where a, (3, 7, k and fx are constants, possess admissible transformations that are 
not generated by transformations from the equivalence group G^ p . The subclass 
of the class ([6]) ; singled out by the condition (fTT|) . is closed under any admissible 
transformation within the class The (constant) parameters of the representa- 
tion (jlip are transformed by an admissible transformation in the following way: 

a = a, (3 = /3, k = 8\/TtK — X if (a, j3) ^ (0, 0), 

7 Ttt ~ V lT tt ~ X t 6X 

li = — — 777=0, A = —jiX — — + 



' T t Tf r T t 2T ( 2 ' r T t y/T t ' 

In particular, Tt = 1/7 ^ 1\x. 

Finally, we can formulate the assertion on group classification with respect to 
the set of admissible transformations. 

Theorem 4. Up to point equivalence, cases of extension of the maximal Lie 
invariance algebras in class © are exhausted by cases 0, 2, 3, 4^ =0 and 6| a=o °S 
tabled 
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3 Algorithm of finding reduction operators 
via mappings between classes 

At first we adduce the definition of nonclassical symmetries [7lll91l26]. adopt- 
ing it for the case of one second-order PDE with two independent variables, 
relevant for this paper. Consider a second-order differential equation C of the 
form L(t, x,U(2)) = for the unknown function u of the two independent vari- 
ables t and x, where um) = (u,ut,u x ,uu,ut x ,u xx ). Let Q be a first-order differ- 
ential operator of the general form 

Q = r(t, x, u)d t + £(t, x, u)d x + T)(t, x, u)d u , (r, £) / (0, 0). 

Definition 1. The differential equation C is called conditionally invariant with 
respect to an operator Q if the relation 

Q (2) L(t,x,u (2) )| £nQ(2) = (12) 

holds, which is called the conditional (or nonclassical) invariance criterion. Then 
Q is called conditional symmetry (or nonclassical symmetry, Q-conditional sym- 
metry or reduction operator) of the equation C. 

The symbol Qr 2 ) stands for the standard second prolongation of Q (see e.g. [HI 
I16j). Q^ 2 ) is the manifold determined in the second-order jet space by the dif- 
ferential consequences of the characteristic equation Q[u] := n — rut — £u x = 0, 
which have, as differential equations, orders not greater than two. 

It was proved in [2B] that a differential equation C is conditionally invariant 
with respect to the operator Q if and only if the ansatz constructed with this 
operator reduces the equation C That is why it seems natural to call operators 
of conditional (nonclassical symmetries) reduction operators. 

Here we present the algorithm of application of equivalence transformations, 
gauging of arbitrary elements and mappings between classes of equations to clas- 
sification of reduction operators of class ([2]) 

1. At first we gauge class ([2]) to subclass © constrained by the condition / = g. 
Then class (J3|) is mapped to the imaged class ([6]) by transformation ([5]). 

2. Reduction operators should be classified up to the equivalence relations gen- 
erated by the corresponding equivalence groups or even by the whole sets of 
admissible transformations. As the singular case r = is "no-go" |10|,l25j . 
only the regular case t / (reduced to the case r = 1) should be considered. 
Operators equivalent to Lie symmetry ones should be neglected. 

3. It is well-known (see e.g. [fl[3]) that the equations from the imaged class © 
with F = and H = const and, therefore, all equations similar to them 
with respect to point transformations possess no regular reduction operators 
that are inequivalent to Lie symmetry operators. This is why all the above 
equations should be excluded from the consideration. 
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4. Preimages of the nonclassical symmetries obtained and of equations admit- 
ting them should be found using the inverses of gauging transformations 
and the push-forwards by these inverses on the sets of operators. 

Reduction operators of equations from class ([3]) are easily found from reduction 
operators of corresponding equations from ([6]) using the formula 

Q = rd t + id x + {r l -^G x )d u . (13) 

Here r, £ and rj respectively are the coefficients of dt, d x and d v in a reduction 
operator of an equation from class 0. The function G is defined in (jHJ). 

In [23 , 24J we discussed two ways to use mappings between classes of equations 
in the investigation of reduction operators and their usage to finding exact so- 
lutions. The preferable way is based on the implementation of reductions in the 
imaged class and preimaging of the obtained exact solutions instead of preimaging 
the corresponding reduction operators. 



4 Reduction operators and exact solutions 

Following the above algorithm, we look for G~ -inequivalent reduction operators 
with nonvanishing coefficient of dt for the equations from the imaged class ([6]). 
Up to the usual equivalence of reduction operators we need to consider only the 
operators of the form 

Q = d t + £(t,x,v)d x + rj(t,x,v)d v . 

Applying conditional invariance criterion (j!2j) to equation (|6|) we obtain a third- 
degree polynomial of v x with coefficients depending on t, x and v, which has to 
identically equal zero. Splitting it with respect to different powers of v x results in 
the following determining equations for the coefficients £ and rj: 

& - ixx + 26^ + 3& {H + ee v ) + 2 Vvx - 2^r? + F£ x + £F X = 0, (14) 
m - Vxx + = i R x + Fr, x + {2i x - r] v )H + ee v (rj + 2£ x - Vv ) . 

Integration of the first two equations of (PT41) gives us the expressions for £ and 
r] with an explicit dependence on v: 

£ = av + 6, 

rj = — a 2 v 3 + (a x — ah — aF)v 2 + cv + d, 
3 

where a = a(t, x), b = b(t, x), c = c(t, x) and d = d(t, x) are smooth functions of t 
and x. 

Substituting the expressions (fl~5j) for £ and rj into the third and forth equations 
of (|14p and collecting the coefficients of different powers of v in the resulting 
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equations, we derive the conditions a = c = 0, d = —2b x and two classifying 
equations, which contain both the coefficient b = b(t, x) and the arbitrary elements 
F = F(x) and H = H(x). Summing up the above consideration, we have the 
following assertion. 

Proposition 2. Any regular reduction operator of an equation from the imaged 
class © is equivalent to an operator of the form 

Q = d t + bd x -2b x d v , (16) 

where the coefficient b = b(t, x) satisfies the overdetermined system of partial 
differential equations 

b t -b xx + 2bb x + Fb x + bF x = 0, 

(17) 

bH x + 2b x H - 466x8. - 2(Fb) xx - 2Fb xx = 

with the corresponding values of the arbitrary elements F = F(x) and H = H(x). 

The second equation of (|17p can be written in the more compact form 

4(6 + F)b xx = 2Kb x + K x b, 

where K = H — 2F X , which is more convenient for the study of compatibility. 

Analogously to the power case, we were not able to completely study all the 
cases of integration of system (|17p depending on values of F and H. This is why 
we try to solve this system under different additional constraints imposed either 
on 6 or on (F, H) . 

The most interesting results are obtained for the constraint 6j = 0. Then F 
and H are expressed, after a partial integration of (fT7|) . via the function 6 = b(x) 
that leads to the following statement. 

Theorem 5. For an arbitrary smooth function b = b(x) the equation from class © 
with the arbitrary elements 

F = l{b x + h- b 2 ) , 

2 (18) 
H = -^(k 2 + b x {h - b 2 ) + bb xx ) , 

where k% and k 2 are constants, admits the reduction operator f)16|) with the same b. 
An ansatz constructed by the reduction operator (|16p with bt = has the form 

v = z(co) — 21n |6|, where to = t — 

The substitution of the ansatz into equation © leads to the reduced ODE 

Zwu ~ fazu + ee z + 2k 2 = 0. (19) 
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For ki = the general solution of (|19p is written in the implicit form 

(ci -Ak 2 z-2ee z )~^dz = ±(w + c 2 ). (20) 



Up to similarity of solutions of equation Q , the constant c 2 is inessential and can 
be set to equal zero by a translation of u, which is always induced by a translation 
of r. 

Setting additionally k 2 = in ([20]) , we are able to integrate in (f20l) in a closed 
form and to explicitly write down the general solution of (|19p . If e = 1 then 
ci > and (|20j) gives the following expression for e z : 

2 s 2 



cosh (siw + s 2 ) 



2s 2 



Here and below si = y |ci|/2 and s 2 = c 2 si. If e = — 1, the integration leads to 

ci > 0, 



sinh {s\UJ + s 2 ) 



2s? 



cos 2 (siw + s 2 ) 
2 



ci < 0, 

a = 0. 



l(^ + c 2 ) 2 ' 

As a result, for the equation from class © of the form 



1 2 

i>t = v xx + - (b x - b 2 ) v x + ee v + - (b x 



bb x ) 



(21) 



with e 



-1 we construct three families of exact solutions 
-2 In 



-2 In 



-2 In 



y/2 , . , ( fdx 
- — 6smh sir — si / — + s 2 
2si V J 

V2 , / fdx 
— 6cos Ut- Sl J T +s 2 



(22) 



V2 



b t 



dx 



+ c 2 



where s\, s 2 and c 2 are arbitrary constants, s% 7^ 0. Also we obtain a family of 
exact solution 



-2 In 



V2 
2si 



b cosh I sir — si 



dx 



+ S 2 



(23) 



of the equation ([2T|) with e = 1. 
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We continue the consideration by studying whether the equations from class ([6]) 
possessing nontrivial Lie symmetry properties, i.e. having the maximal Lie invari- 
ance algebras of dimension two or three, have nontrivial (i.e., inequivalent to Lie 
ones) regular reduction operators. It has been already remarked that constant 
coefficient equations from class © do not admit such reduction operators [HE]. 
Hence it is needless to consider cases 4 and 6 of table 1 as well as case 5 connected 
with case 6 by point transformation (jlOp . As case 1 reduces to case 2 with the 
same transformation (jlOH . we have to study only two cases, namely cases 2 and 3. 
We substitute the pairs of values of the parameter-functions F and H correspond- 
ing to cases 2 and 3 into system (|17p in order to find relevant values for b. We 
ascertain that bt = is a necessary condition for existing non-Lie regular reduc- 
tion operators for equations with the above values of (F, H). This is why we can 
use equations (fTHj) instead of (fTTl) for further studying. 

The investigation of case 3 of table 1 leads to the conclusion that there are no 
non-Lie regular reduction operators for this case. 

The functions F and H presented in case 2 of table 1 satisfy (118p if and only 
if f3 = 2(1 — a), i.e. they have the form F = ax^ 1 , H = 2(1 — a)x~ 2 , and 
ki = &2 = 0. The corresponding value of b is b = —(1 + a)x _1 . Hence a ^ — 1 
since otherwise b = 0. Substituting the derived form of the function b into the 
formulas (i22"j) and (f23|) . we find that the equation 



a v 2(1 -a) 
v t = v xx + -v x + ee H 3 

X X A 

has the families of exact solutions 



(24) 



-2 In 



V2(l + 



a) 



2s\x 



cosh sit + 



S\X 



2(1 + a] 



+ S 2 



if e = 1 and 



-2 In 



-2 In 



-2 In 



2s\x 
y/2(l + a 



2s\x 
V2(l + a) 



2x 



COS Sit + 



t + 



2(1 + a) 



2(1 + a) 

2 

SlX 

2(1 + a) 

+ C 2 



+ S 2 



+ S 2 



if e = —1. Recall that s\, S2 and C2 are arbitrary constants, si + 1 0. 

As a representative of the preimage of equation (|24p with respect to the trans- 
formation ([5]), we can choose the equation 



x a u t = (x a u x ) x + ex a+2 e u . 



(25) 
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Solutions of this equation can be easily constructed from the above solutions of 
equation (I24p using the transformation u = v — 2 hi \x\. If a = 1, the chosen 
equation (j25|) can be replaced, e.g., by xu t = {xu x ) x + exe u which is just another 
representation of equation (|24|) . 

Non-Lie exact solutions of the equation 



where a / -1 (case 1 of table 1), can be easily obtained from exact solutions 
of the equation (j25j) using the transformation (jlOp . The corresponding reduction 
operator has the form (|16H with b = —(1 + a)x _1 — fix. 
We also prove the following assertions. 

Proposition 3. Equations from class ([6]) with F = const or H = const may 
admit only nontrivial regular reduction operators that are equivalent to operators 
of the form (j!6j) . where the function b does not depend on the variable t. 

Proposition 4. Any reduction operator of an equation from class ([6]), having the 
form (|16p with b xx = 0, is equivalent to a Lie symmetry operator of this equation. 

The proofs of these propositions are quite cumbersome and will be presented 
elsewhere. 
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